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The observed late-time acceleration of the cosmic expansion constitutes a fundamental problem 
in modern theoretical physics and cosmology. In an attempt to weight the validity of a large number 
of dark energy models, I use the recent measurements of the expansion rate of the Universe, the 
clustering of galaxies the CMB fluctuations as well as the large scale structure formation, to put 
tight constraints on the different models. 
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1. INTRODUCTION 

Recent studies in observational cosmology lead to the 
conclusion that the available high quality cosmological 
data (Type la supernovae, CMB, etc.) are well fitted by 
an emerging "standard model" . This standard model, 
assuming flatness, is described by the Friedman equation: 



H 2 (a) 



IttG 



[p m {a) + Px(a)} 



(1) 



where a(t) is the scale factor of the universe, p m (a) is the 
density corresponding to the sum of baryonic and cold 
dark matter, with the latter needed to explain clustering, 
and an extra component px{a) with negative pressure, 
called dark energy, needed to explain the observed accel- 
erated cosmic expansion (eg., Davis et al. 2007; Kowalski 
et al. 2008; Komatsu et al. 2009 and references therein). 

The nature of the dark energy is one of the most 
fundamental and difficult problems in physics and cos- 
mology. Indeed, during the last decade there has been 
a theoretical debate among the cosmologists regarding 
the nature of the exotic "dark energy" . Various candi- 
dates have been proposed in the literature, such as a cos- 
mological constant A (vacuum), time- varying A(t) cos- 
mologies, quintessence, k— essence, vector fields, phan- 
tom, tachyons, Chaplygin gas and the list goes on (for 
recent reviews see Peebles & Ratra 2003; Copeland, Sami 
& Tsujikawa 2006; Caldwell & Kamionkowski 2009). 
Within this framework, high energy field theories gener- 
ically indicate that the equation of state of such a dark 
energy is a function of cosmic time. To identify this type 
of evolution of the equation of state, a detailed form of 
the observed H(a) is required which may be obtained 
by a combination of multiple dark energy observational 
probes. 

In this cosmological framework, a serious issue is how 
(and when) the large scale structures form. Galaxies 
and large-scale structure grew gravitationally from tiny, 
nearly scale-invariant adiabatic Gaussian fluctuations. In 
this paper we focus also on galaxy clusters that occupy 
an eminent position in the structure hierarchy, being the 
most massive virialized systems known and therefore they 
appear to be ideal tools for testing theories of struc- 
ture formation and extracting cosmological information. 



The cluster distribution basically traces scales that have 
not yet undergone the non-linear phase of gravitation- 
ally clustering, thus simplifying their connections to the 
initial conditions of cosmic structure formation. 

The structure of the paper is as follows. In section 2 
we briefly discuss the dark energy issue. In sections 3 and 
4 we present the various dark energy models and we use 
a joint statistical analysis, in order to place constraints 
on the main cosmological parameters. In section 5 we 
present the corresponding theoretical predictions regard- 
ing the formation of the galaxy clusters. Finally, we draw 
our conclusions in section 6. Throughout the paper we 
will use H ~ 71km/sec/Mpc. 



2. THE DARK ENERGY EQUATION OF STATE 

In the context of general relativity it is well known that 
for homogeneous and isotropic fiat cosmologies, driven by 
non-relativistic matter and dark energy with equation of 
state px — w(a)px, the expansion rate of the Universe 
can be written as (see eqU]) 



E 2 (a) 



H 2 (a) 
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Note, that E(a) is the normalized Hubble flow, fi m is 
the dimensionless matter density at the present epoch, 
fix = 1 — Q m is the corresponding dark energy density 
parameter and w(a) is the dark energy equation of state. 
Inverting the above equation we simply derive 
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The exact nature of the dark energy has yet to be found 
and thus the dark energy equation of state parameter 
includes our ignorance regarding the physical mechanism 
which leads to a late cosmic acceleration. 

On the other hand, it is possible to extent the previous 
methodology in the framework of modified gravity (see 
Linder & Jenkins 2003; Linder 2004). In this scenario, it 
is assumed that the dark energy may be an illusion, indi- 
cating the need to revise the general relativity and thus 
also the Friedmann equation. From the mathematical 
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point of view, it can be shown that instead of using the 
exact Hubble flow through a modification of the Fried- 
mann equation we can utilize a Hubble flow that looks 
like the nominal one (see eqUJ. The key point here is 
to consider that the accelerated expansion of the uni- 
verse can be attributed to a "geometrical" dark energy 
component. Due to the fact that the matter density in 
the universe (baryonic+dark) can not accelerate the cos- 
mic expansion, we perform the following parametrization 
(Linder & Jenkins 2003; Linder 2004): 



E\a) 



H 2 (a) 
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(4) 



Obviously, with the aid of the latter approach we include 
any modification to the Friedmann equation of general 
relativity in the last term of eq.(H]). Now from eqs.([31 
H]) we can derive, after some algebra, the "geometrical" 
dark energy equation of state 



_ 1 f dlnSH 2 
3 eflna 



(5) 



From now on, for the modified cosmological models we 
will use the above formulation. 



from the WMAP 5-years data (Komatsu et al. 2009) is 
R = 1.71 ± 0.019 at z u = 1090 [or a ts = (1 + z/s)" 1 ~ 
9.17x 10~ 4 ]. In this case, the Xcmb function is given 



Xcmb(P) 



[R(p) 1-71] ; 
0.019 2 



(7) 



Finally, we use the Union08 sample of 307 supernovae 
of Kowalski et al. (2008). The corresponding XsNia func- 
tion becomes: 



XSNIa(P) 



i=l 



(8) 



where a, = (1 + Zi)^ 1 is the observed scale factor of the 
Universe, is the observed redshift, is the distance 
modulus fj, = m — M = 51og<ii + 25 and di(a,p) is 

the luminosity distance <ii(a,p) = ^ J a y 2H( y ) where c 
is the speed of light. We can combine the above probes 
by using a joint likelihood analysis: C to t(p) = £bao x 

£cmb X £ S NIa Or Xfot(P) = X.BAO + Xcmb + XSNIa m 0rder 

to put even further constraints on the parameter space 
used. Note, that we define the likelihood estimator 1 as: 
Cj oc exp[-x|/2]. 



3. LIKELIHOOD ANALYSIS 

We will use various cosmological observations in order 
to constrain the dark energy models described in section 
4. First of all, we use the Baryonic Acoustic Oscilla- 
tions (BAOs). BAOs are produced by pressure (acoustic) 
waves in the photon-baryon plasma in the early universe, 
generated by dark matter overdensities. Evidence of this 
excess was recently found in the clustering properties of 
the luminous SDSS red-galaxies (Eisenstein et al. 2005; 
Padmanabhan et al. 2007) and it can provide a "stan- 
dard ruler" with which we can constraint the dark energy 
models. In particular, we use the following estimator: 

- 2/3 

A(p) = 



measured from the 



[z2£(a s )]i/3 

SDSS data to be A = 0.469 ± 0.017, where z s = 0.35 [or 
a s = (1 + Zg)^ 1 ~ 0.75]. Therefore, the corresponding 
Xbao function is simply written 



Xbao(p) 



[A(p) - 0.469] 2 
0.017 2 



(6) 



where p is a vector containing the cosmological parame- 
ters that we want to fit. 

On the other hand, a very accurate and deep geo- 
metrical probe of dark energy is the angular scale of 
the sound horizon at the last scattering surface as en- 
coded in the location lf T of the first peak of the Cos- 
mic Microwave Background (CMB) temperature per- 
turbation spectrum. This probe is described by the 
CMB shift parameter (Bond, Efstathiou & Tegmark, 
1997; Nesseris & Perivolaropoulos 2007) and is defined 
as R = V^m f ai a ' id E{a) - ^ ne sn ^ parameter measured 



4. CONSTRAINTS ON THE FLAT DARK 
ENERGY MODELS 

In this section, we consider a large family of flat dark 
energy models and with the aid of the above cosmologi- 
cally relevant observational data, we attempt to put tight 
constraints on their free parameters. In the following 
subsections, we briefly present these cosmological mod- 
els which trace differently the nature of the dark energy. 



4.1. Constant equation of state - QP model 

In this case the equation of state is constant (see for 
a review, Peebles & Ratra 2003; hereafter QP-models). 
Such dark energy models do not have much physical mo- 
tivation. In particular, a constant equation of state pa- 
rameter requires a fine tuning of the potential and ki- 
netic energies of the real scalar field. Despite the lat- 
ter problem, these dark energy models have been used 
in the literature due to their simplicity. Notice that 
dark energy models with a canonical kinetic term have 
— 1 < w < —1/3. On the other hand, models of phantom 
dark energy (w < —1) require exotic nature, such as a 



1 Likelihoods are normalized to their maximum values. Note, that 
the step of sampling is 0.01 and the errors of the fitted parameters 
represent 2a uncertainties. Note that the overall number of data 
points used is Ntot = 309 and the degrees of freedom: dof = 
Ntot — "fit, with ngt the number of fitted parameters, which 
vary for the different models. 
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scalar field with a negative kinetic energy. Now using 
eq.© the normalized Hubble parameter becomes 



E 2 (a) = tt m a~ 3 + (1 - tt m )a 



-3(l+«>) 



(9) 



Comparing the QP-models with the observational data 
(we sample Q m G [0.1,1] and w € [—2,-0.4]) we find 
that the best fit values are f2 TO = 0.28 ± 0.02 and 
w = -1.02 ±0.06 with Xtot(^m,w)/dof ~ 309.2/307 
in very good agreement with the 5 years WMAP data 
Komatsu et al. (2009). Also Davis et al. (2007) utilizing 
the Essence-SNIa+BAO+CMB and a Bayesian statis- 
tics found Q m = 0.27 ± 0.04, while Kowalski et al. 
(2008) using the Union08-SNIa+BAO+CMB obtained 



i 1 1 ' 



0.285 



+0.03 
-0.03- 



Obviously, our results coincide within 



la errors. It is worth noting that the concordance A- 
cosmology can be described by QP models with w strictly 
equal to -1. In this case we find: il m = 0.28 ± 0.02 with 
xL(O m )/do/~ 309.3/308. 



4.2. The Braneworld Gravity - BRG model 

In the context of a braneworld cosmology (hereafter 
BRG) the accelerated expansion of the universe can be 
explained by a modification of gravity in which gravity 
itself becomes weak at very large distances (close to the 
Hubble scale) due to the fact that our four dimensional 
brane survives into an extra dimensional manifold (Def- 
fayet, Dvali, & Cabadadze 2002). The interesting point 
in this scenario is that the corresponding functional form 
of the normalized Hubble flow, eq. @, is affected only 
by one free parameter (O m ). Notice, that the quantity 
5H 2 is given by 



m 2 



bw 



(10) 



where Clbw = (1 — ^?n) 2 /4- The geometrical dark energy 
equation of state parameter (see eq(5]) reduces to 



1 



wla) = — 



1 + £l m (a) 



(11) 



where J7 m (a) = fl m a~ 3 / E 2 (a) . The joint likelihood anal- 
ysis provides a best fit value of Q m = 0.24 ± 0.02, but the 
fit is rather poor xiot(^m)/dof ~ 369/308. 



4.3. The parametric Dark Energy model 
model 



CPL 



In this model we use the Chevalier-Polarski-Linder 
(Chevallier & Polarski 2001; Linder 2003, hereafter CPL) 
parametrization. The dark energy equation of state 
parameter is defined as a first order Taylor expansion 
around the present epoch: 



w(a) — wq + w%(l — a) 



(12) 



The normalized Hubble parameter is given by (see eqj2]) : 
E 2 (a) = ^ m a- 3 + (l-a»)a~ 3(1+ "' 0+,1 ' l) e 3 " Jl(a " 1) , (13) 

where Wq and u>i are constants. We sample the un- 
known parameters as follows: wq G [—2, —0.4] and w± G 
[-2.6,2.6]. We find that for the prior of Q m = 0.28 the 
overall likelihood function peaks at wq — — l.ll ; 16 and 
W\ = 0.60^5'54 while the corresponding Xtot( w o, w\)/dof 
is 307.6/307^ 



4.4. The low Ricci dark energy - LRDE model 

In this modified cosmological model, we use a sim- 
ple parametrization for the Ricci scalar which is based 
on a Taylor expansion around the present time: 1Z = 
fo + — a) [for more details see Linder 2004]. It is 
interesting to point that at the early epochs the cosmic 
evolution tends asymptotically to be matter dominated. 
In this framework, we have 



E 2 (a) = 



a 4(r +ri-l) e 4ri(l-o) 



a> a t 
a < a t 



(14) 



where a t = 1 — (1 — 4Yo)/4ri. The matter density pa- 
rameter at the present epoch, is related with the above 



constants via fi m = 



_ / 4r -4ri-l \ 
\ 4r i J 



4r +4ri — 1 



The 



equation of state parameter that corresponds to the cur- 
rent geometrical dark energy model is given by 



w(a) = 



I -ATI 



(15) 



Note, that we sample the unknown parameters as follows: 
r G [0.5,1.5] and r\ G [—2.4,-0.1] and here are the 



results: 



ro 



0.82 ± 0.04 and r x = -0.74 



+0.10 
-0.08 



0.28) with Xtot( r o^x)/dof ~ 309.8/307. 



4.5. The high Ricci dark energy - HRDE model 

A different than the previously described geometrical 
method was defined by Linder & Cahn (2007), in which 
the Ricci scalar at high redshifts evolves as 



1 — 3wi ■ 



SH 2 
H 2 



(16) 



where SH 2 = E 2 (a)-Q r , 



,-3 



In this geometrical pattern 



the normalized Hubble flow becomes: 

E 2 (a) = Q m a- 3 (i + (Sa -3» 1 )- 1 ^/Mi+« ^ 

where /? = fl^ 1 — 1. Using the same sampling as in the 
QP-modcls, the joint likelihood peaks at f2 m = 0.28 ± 
0.03 and w x = -1.02 ± 0.1 with xLt(^ m ,wi)/dof ~ 
309.2/307. To this end, the effective equation of state 
parameter is related to E(a) according to eq.©. 
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4.6. The tension of cosmological magnetic fields 
TCM model 



Recently, Contopoulos & Basilakos (2007) proposed a 
novel idea which is based on the following consideration 
(hereafter TCM): if the cosmic magnetic field is gener- 
ated in sources (such as galaxy clusters) whose overall 
dimensions remain unchanged during the expansion of 
the Universe, the stretching of this field by the cosmic 
expansion generates a tension (negative pressure) that 
could possibly explain a small fraction of the dark en- 
ergy (~ 5 — 10%). In this fiat cosmological scenario the 
normalized Hubble flow becomes: 



E\a) 



-3+n 



(18) 



where fis is the density parameter for the cosmic mag- 
netic fields and fiA = 1 — fim ~ The equation of 
state parameter which is related to magnetic tension is 
(see eq[3]) 



w(a) 



3fi/ 



-3+n 



(19) 



To this end, we sample fi_B S [0,0.3] and n £ [0,10] 
and we find that for the prior of fi m = 0.28 the best 
fit values are: fi B = 0.10 ± 0.10 and n = 3. 60^1 with 
xL(O s ,n)/do/~ 308.9/307. 



4.7. The Pseudo-Nambu Goldstone boson - PNGB 
model 

In this cosmological model the dark energy equation 
of state parameter is expressed with the aid the poten- 
tial V((j>) oc [1 + cos(<j)/F)} (Abrahamse et al 2008 and 
references therein): 



w(a) =—! + (! + wo)a F , 



(20) 



where w G [-2,-0.4] and F € [0,8]. In case of a <C 1 
we get w(a) — > —1. Based on this parametrization the 
normalized Hubble function is given by 

E 2 {a) = fi m a" 3 + (1 - n m ) Px {a) . (21) 

In this context, the corresponding dark energy density is 

"3(1 + w ), 



px{a) = exp 



-il-a*) 



(22) 



Notice, that the likelihood function peaks at wq 
-1.04 ±0.17 and F = 5.9 ± 3.2 with Xtot( w 0, F)/dof 
317/307. 



4.8. The early dark energy - EDE model 



The basic assumption here is that at early epochs the 
amount of dark energy is not negligible (Doran, Stern & 
Thommes 2006 and references therein). In this model, 
the total dark energy component is given by: 



n x ( a ) 



1 - fir, 



n e (i 



-3w 



1 fim fir, 



fie(l 



-3wo 



), 

(23) 

where fi e is the early dark energy density and Wq is the 
equation of state parameter at the present epoch. Notice, 
that the EDE model was designed to simultaneously (a) 
mimic the effects of the late dark energy and (b) provide 
a decelerated expansion of the universe at high redshifts. 
The normalized Hubble parameter is written as: 



E\a) 



\tqrnCL 



l-fix(o) ' 



(24) 



while using eq.Q, we can obtain the equation of state 
parameter as a function of the scale factor. From the 
joint likelihood analysis we find that fi e = 0.05 ± 0.04 
and w Q = -l.Utoll ( for thc P rior of fi m = °- 28 ) with 
x 2 tot (n e ,w )/dof ~ 308.7/307. 



4.9. The Variable Chaplygin Gas - VCG model 

Let us consider now a completely different model 
namely the variable Chaplygin gas (hereafter VCG) 
which corresponds to a Born-Infeld tachyon action 
(Bento, Bertolami & Sen 2004). Recently, an interesting 
family of Chaplygin gas models was found to be consis- 
tent with the current observational data (Dev, Alcaniz & 
Jain 2003). In the framework of a spatially flat geome- 
try, it can be shown that the normalized Hubble function 
takes the following formula: 



E 2 (a) = n b a- 3 + (l-n b )^/B s 



(l-B s )a-™ , (25) 



where fif, ~ 0.021/i~ 2 is the density parameter for the 
baryonic matter (see Kirkman et al. 2003) and B s E 
[0.01,0.51] and n e [-4,4]. The effective matter density 
parameter is: fi^ = fi&+(l — fi&)\/B7. We find that the 
best fit parameters are B s = 0.05±0.02 and n = 1.58to43 
(Sl%f ~ 0.26) with Xtot(B s ,n)/dof ~ 314.7/307. 



5. EVOLUTION OF MATTER 
PERTURBATIONS 

Thc evolution equation of the growth factor for models 
where the dark energy fluid has a vanishing anisotropic 
stress and the matter fluid is not coupled to other matter 
species is given by (Linder & Jenkins 2003) 



Another cosmological scenario that we include in our 
paper is the early dark energy model (hereafter EDE). 



D 



1 - 



w(a) 
1 + X(a) 



D 



3 X(a) D 
21 + X(a) a? 



= , (26) 



5 



where 

X(a) = -^-e- 3 ^ 1 = . (27) 

Note, that the prime denotes derivatives with respect to 
the scale factor. Useful expressions of the growth factor 
can be found for the ACDM cosmology in Peebles (1993), 
for dark energy models with w = const in Silveira, & 
Waga (1994) for dark energy models with a time vary- 
ing equation of state in Linder & Cahn (2007) and for 
the scalar tensor models in Gannouji & Polarski (2008). 
In this work the growth factor evolution for the current 
cosmological model is derived by solving numerically eq. 
(|2"oT) . Note, that the growth factors are normalized to 
unity at the present time. 



5.1. The formation of galaxy clusters 

In this section we briefly investigate the cluster forma- 
tion processes by generalizing the basic equations which 
govern the behavior of the matter perturbations within 
the framework of the current dark energy models. Also 
we compare our predictions with those found for the tra- 
ditional A cosmology. This can help us understand better 
the theoretical expectations of the dark energy models as 
well as the variants from the usual A cosmology. 
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FIG. 1: The predicted fractional rate of cluster formation 
as a function of redshift for the current cosmological models 
(ag — 0.80). The points represent the following cosmological 
models: (a) BRG (open stars), (b) LRDE (open squares), (c) 
TCM (open triangles), (d) EDE (open circles) and (e) PNGB 
(solid triangles). The lines represent: (a) CPL model (long 
dashed) , (b) HRDE model (dot line) and VCG model (dashed 
line) . 

The concept of estimating the fractional rate of clus- 
ter formation has been proposed by different authors (eg., 
Weinberg 1987; Richstone, Loeb & Turner 1992). In par- 
ticular, these authors introduced a methodology which 
computes the rate at which mass joins virialized struc- 
tures, which grow from small initial perturbations in the 
universe. The basic tool is the Press & Schechter (1974) 
formalism which considers the fraction of mass in the 



universe contained in gravitationally bound structures 
(such as galaxy clusters) with matter fluctuations greater 
than a critical value S c , which is the linearly extrapolated 
density threshold above which structures collapse (Eke, 
Cole & Frenk 1996). Assuming that the density con- 
trast is normally distributed with zero mean and variance 
a 2 (M, z) we have: 



V(5,z) 



2tt<t(M, z 



■exp 



2cr 2 (M,; 



(28) 



In this kind of studies it is common to parametrize the 
rms mass fluctuation amplitude at 8 h Mpc which can 
be expressed as a function of redshift as cr(M, z) — 
crg(z) — D(z)ag. The current cosmological models are 
normalized by the analysis of the WMAP 5 years data 
erg = 0.80 (Komatsu et al. 2009). The integration of 
eq. ([28l) provides the fraction of the universe, on some 
specific mass scale, that has already collapsed producing 
cosmic structures (galaxy clusters) at redshift z and is 
given by Richstone et al. (1992): 



P(z) 



V(6,z)dS=- 



1-erf 



V2a 8 (z) 



(29) 

Notice, that for the model of modified gravity (BRG) 
we use S c ~ 1.47 (see Schafer & Koyama 2008), for the 
EDE model we use S c ~ 1.4 (see Bartelmann, Doran & 
Wetterich 2006). For the rest of the dark energy models, 
due to the fact that w ~ — 1 close to the present epoch, 
we utilize a S c approximation which is given by Weinberg 
& Kamionkowski (2003 see their eq.18). 

Obviously the above generic form of ea. ([29)) depends 
on the choice of the background cosmology. The next 
step is to normalize the probability to give the number 
of clusters which have already collapsed by the epoch z 
(cumulative distribution), divided by the number of clus- 
ters which have collapsed at the present epoch (z = 0), 
F(z) = P(z) / P(0). In figure 1 we present, in a logarith- 
mic scale, the behavior of normalized cluster formation 
rate as a function of redshift for the various dark energy 
models. In general, prior to z ~ 0.2 the cluster formation 
has terminated due to the fact that the matter fluctua- 
tion field, D(z), effectively freezes. For the traditional A 
cosmology we find the known behavior in which galaxy 
clusters appear to be formed at high redshifts z ~ 2 
(Basilakos 2003; Basilakos, Sanchez & Perivolaropoulos 
2009 and references therein). From figure 1 it becomes 
also clear, that the vast majority of the dark energy mod- 
els seem to have a cluster formation rate which is close to 
that predicted by the usual A cosmology (see solid line 
in figure 1). However, for the BRG and EDE cosmologi- 
cal scenarios we find that galaxy clusters appear to form 
earlier (z ~ 3) than in the CPL, LRDE, HRDE, TCM, 
PNGB and VCG dark energy models. 
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6. CONCLUSIONS 

In this work we have studied analytically and numer- 
ically the overall dynamics of the universe for a large 
number of dark energy models beyond the concrdance 
A cosmology, by using several parameterizations for the 
dark energy equation of state. We performed a joint like- 
lihood analysis, using the current high-quality observa- 



tional data (SNIa, CMB shift parameter and BAOs), and 
we put tight constraints on the main cosmological param- 
eters. We also find that for the vast majority of the dark 
energy models, the large scale structures (such as galaxy 
clusters) start to form prior to z ~ 1 — 2. 
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